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Abstract 


Recently Csaki, Kaloper and Terning (|hep-ph/0111311|) snggested that the ob¬ 
served dimming of distant type la supernovae may be a consequence of mixing of the 
photons with very light axions. We point out that the effect of the plasma, in which 
the photons are propagating, must be taken into account. This effect changes the 
oscillation probability and renders the dimming frequency-dependent, contrary to 
observations. One may hope to accommodate the data by averaging the oscillations 
over many different coherence domains. We estimate the effect of coherence loss, 
either due to the inhomogeneities of the magnetic field or of the intergalactic plasma. 
These estimates indicate that the achromaticity problem can be resolved only with 
very specihc, and probably unrealistic, properties of the intergalactic medium. 


1 Introduction 


The mixing between photons and axions, in an external magnetic held, is a well 
studied mechanism i I, I, i as is its analog with gravitons p, |^. It is exper¬ 
imentally used, since the pioneer work by Sikivie [0, to constrain the axion pa¬ 
rameters ||7|, PI, |9|, pT, 0 (see also e.g. Ref. |]^ for an up to date review on such 
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experiments). The astrophysical and cosmological implications of this mechanism 
have also been studied 13. and it was recently advocated to be a possible expla¬ 
nation for the observed dimming of distant type la supernovae [p!3| , by Csaki et 
B e also [0 for the analysis of supernovae data with the model of |T5[|). 


al 

The underlying idea is simply that the luminosity of distant supernovae can be di¬ 
minished due to the decay of photons into very light pseudoscalar particles induced 
by intergalactic magnetic helds over cosmological distances, and thus that the lumi¬ 
nosity distance-redshift relationship can mimic the one of a universe with a non zero 
cosmological constant without need for a cosmological constant. The pseudoscalar 
particles must have an electromagnetic coupling similar to axions, and a specihc and 
very small mass (m ~ 10“^® eV) to avoid affecting the cosmic microwave background 
anisotropy beyond its observed value. Another aspect of photon-pseudoscalar con¬ 
version in intergalactic magnetic helds is the change of the polarization properties 
of distant sources 0 such as supernovae. 

The implications on the cosmic microwave background of the similar effect in¬ 
volving photon-graviton oscillation has been considered in empty space 10, 0] and 
it was shown that it becomes negligible for standard cosmological magnetic helds 
once the contribution of the intergalactic plasma is properly taken into account (the 
case of axions is also considered in |^). The ehect of the inhomogeneities of the 
electron density upon the coherence of the oscillations was also considered by Carl¬ 
son et al. 0 . 

The ehect of this plasma was not addressed in [0, and it is our purpose to 
discuss it in this work. 

The paper is organized as follows: we hrst remind standard results on photon- 
pseudoscalar oscillations in order to introduce our conventions (see |l|], or 0| where 
contributions from Kaluza-Klein modes were also included). We then discuss specif¬ 
ically the ehect of the plasma for the parameters relevant for type la supernovae. 


2 Photon—pseudoscalar mixing 


We consider the generic action for the pseudoscalar-photon system 


^4 = 


I d*i [-i {S'‘aS„a + my} + 


( 1 ) 


where a is the pseudoscalar held and ttIq its mass. = e^i.po-T’^'^/2 is the dual 
of the electromagnetic tensor, being the completely antisymmetric tensor such 
that 60123 = +1- The pseudoscalar couples to the photon with the coupling IjMa- 
We now consider an electromagnetic plane wave in the presence of a magnetic 
held Bq which is assumed constant on a characteristic scale Ac in the sense that 
its variations in space and time are negligible on scales comparable to the photon 
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wavelength and period. From the three independent vectors 



5, 


O.n 




e±, 


0,n 


(2) 


we dehne a direct orthonormal basis of the three dimensional space {e\\, e^, n 

is the component of Bq perpendicular to the direction of propagation k. 

The electromagnetic wave derives from a potential vector that can be chosen 
to be of the form A = z(74||(s), y4_L(s), 0)e“*“* where s is the coordinate along the 
direction of propagation^. With this decomposition, the coupled Klein-Gordon and 
Maxwell equations derived from the action (|^) read 

DA = (3) 


where A =||, T denotes the polarization. 

Assuming that the magnetic held varies in space on scales much larger than the 
photon wavelength, we can perform the expansion o;^ + = (a; + ids){uj — ids) = 

[uj + k){uj — ids) for a held propagating in the +s direction. If we assume a general 
dispersion relation of the form uj = nk and that the refractive index n satishes 
|n —1| -C 1, we may approximate uj + k = 2uj and k/uj = 1. This approximation can 
be understood as a WKB limit where we set A(s) = |A(s)|e®^^ and assumes that the 
amplitude |A| varies slowly, i.e. that 9^1 A| -C k\A\. In that limit, the above system 
(I) reduces to the linearized system 


(a; — ids + M.) 


All 

A|| 

a 


= 0 . 


The mixing matrix M. is dehned by 


M = 


Ai 0 

0 M\ 


with 


-Mil = 


A|| A 
Am A^ 


M 


The coefficients Am and A^ are given by 


A 2i?o,n 

Am — 




™ 2ci' 


( 4 ) 


( 5 ) 


( 6 ) 


The terms Aa can be decomposed as Aa = Aqed + Aqm + Apiasma- The hrst term 
contains the ehect of vacuum polarization giving a refractive index to the photon 

®with this convention the electric field of the wave is simply E = —dtA = 
M||(s),L^Ti(s),0)e-*“*. 
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(see e.g. Ref. [^) and can be compnted by adding the Enler-Heisenberg effective 
Lagrangian which is the lowest order term of the non-linearity of the Maxwell eqna- 
tions in vacnnm (see e.g. I 


) to the action (|^|j. The second term describes 
the Cotton-Monton effect, i.e. the birefringence of gases and liqnids in presence of 
a magnetic held and the third term the effect of the plasma (since, in general, the 
photon does not propagate in vacnnm). Their explicit expressions are given by 


^QED — 


= 

2 


'^QED ~ 2a;^, 


A 


plasma 




CM 

2 


‘^plasma 

2^’ 

= 2ttCBI 


(7) 

with ^ = (a/457r)(i?o,n/-Bc)^5 = ml/e = 4.41 x 10^^ G, mg the electron mass, 

e the electron charge and a the hne strnctnre constant. C is the Cotton-Monton 
constant ; its effect is to give only the difference of the refractive indices and the 
exact valne of C is hard to determine The plasma freqnency cnpiasma is dehned 
by 


^plasma — 4 : 710 ! 


Uf. 




( 8 ) 


He being the electron density. Note that is always negative whereas Aa is positive 
if the contribntion of the vacnnm dominates and negative when the plasma term 
dominates. 

As seen from Eq. (^, only the component ||, i.e. parallel to the magnetic held, 
conples to the psendoscalars, a hrst conseqnence of which is that the polarization 
plane of a light beam traveling in a magnetic held will rotate. 

The solntion to the eqnation of motion (^ is obtained by diagonalizing Ady 
throngh a rotation 




with the mixing angle'd given by 


cos'd 
— sin'd 


sin'd 

cos'd 


All 


tan 2'd = 2 


A 


M 


All - A, 


(9) 


( 10 ) 


By solving Eq. (^) in this new basis, one can easily compnte the probability of 
oscillation of a photon after a distance of hight s starting from the initial state 
(A||(0) = l,a(0) = 0). It is explicitly given by 

F (7 ^ a) =1 (A||(0) I a(s)) 1^ = sin^ (2^9) sin^ , (11) 

sin^(Aoscs/2) 


= i^M-sY 


( 12 ) 


_ (AoscS/2)2 

®The equation of motion derived from (0) is (^ with A A = 0. We intentionally omit the Euler- 
Heisenberg contribution in the presentation for the sake of clarity. Its Lagrangian is explicitly 
given by Ceh = 


90m* 
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(13) 


with the (reduced) oscillation wavenumber Aqsc given by 

^ _ ^11 ~ _ 2Am 

COS 2'd sin 2-^ 

The oscillation length is thus given by lose = 27r/Aosc- We see that a complete 
transition between a photon and a pseudoscalar is only possible when the mixing is 
maximal (strong mixing regime) i.e. when ■d ~ 7r/4. 


3 Application to Supernovae 


The quantities required for our discussion are Am, A^, Apiasma and Aqed respec¬ 
tively given by equations (H) and (0). It is useful to rewrite them as 


Am 

1 cm“^ 


1 cm ^ 

Apiasma 

1 cm“^ 
Aqed 
1 cm“i 


2 X 






10-9 G; VlQiiGeV 


-1 


-2.5 X 10“2® 
-3.6 X 10-2^ 
1.33 X 10“^® I 


" ma y / UJ \ ^ 

ao-i6evy Uevy ’ 

/ Ue \ 

AeV/ V10-'^cm“3/’ 

0^ \ / Bp 

ievyuo-9Gy ’ 


(14) 


where we have used the facts that leV ~ 5 x 10"^ cm 1 G ~ 1.95 x 10 ^eV^ in 
the natural Lorentz-Heaviside units where a = e^/dvr = 1/137. 


The parameters chosen in are Ma ~ 4 x lO^^GeV, ma ~ 10-^®eV and 
Bq ~ lO-^G. The intergalactic medium (IGM) today is fully ionized, as indicated 
by the lack of Gunn-Petterson effect [^. Thus the mean electronic density can be 
estimated to be (see e.g. 


n. 


~ 10 


-7 


cm 


-3 


(15) 


One immediately sees that, with this choice of parameters, Aqed is always negligible, 
whereas one has | Apiasma | S> |Am| so that the plasma effects are always dominant 
over the pure mass term of the pseudoscalar. 


In order to be more specihe, let us compare the mixing angles d, oscillation 
wavenumbers Aqsc and oscillation probabilities P{'y ^ a) with and without including 
the effect of the intergalactic plasma, and for similar choice of parameters (we will 
use a subscript o for the values computed without the plasma effect). 

When one sets A^ to zero in Eq. (^, as implicitly done in [^, the mixing angle 
reduces to 

t^24--2^~40(^). (16) 
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so that for optical photons (with uj ~ 1.5 eV to 3 eV) tan2'i9o ^ 1 and "do ~ 7r/4. 
This corresponds to a regime in which uj S> m?/A m and the oscillation proba¬ 
bility does not depend on uj so that the oscillation is achromatic. The oscillation 
wavennmber Aqsc.o is given then by 

Aosc.o ~ 2AAr, (17) 


which is also independent of ui. With the choice of parameters nsed in [O, Aqsc.o — 
10 “^®cm“^, so that the oscillation length is larger than the size s of the domain of 
coherence of the magnetic held considered which is of order of a Mpc (~ 3 x1024cm). 
The probability of oscillation over a domain of size s is then well approximated by 


Poil ^ a) {AmsY, (18) 

which is of order 10“4. The nnmber of snch domains in onr Hnbble radius and 
on a given line of sight is given by H^^/s. If one considers that the universe is made 
by patching together such domains with uncorrelated Bq, the coherence is lost from 
domain to domain and one can simply sum up the probability of conversion over 
each domain to obtain the probability of conversion of a photon on cosmological 
distances given by (see also section 4) 

-Po,tot(7 ^ a) ~ AmsHq^. (19) 


This number is of order 1, and one can thus expect a signihcant reduction of the 
luminous hux over cosmological distances. This is the bottom line of the mechanism 
proposed by Csaki et al. [|^ . 

Let us now include plasma effects. Since |Apiasma| S> |Am|, the mixing angle -d 
is now much smaller than -do and 




A 


M 


A 


plasma 


( 20 ) 


of order 10 ^ — 10 2 . With such a low mixing angle {weak mixing regime), the 
probability of oscillation P{'y ^ a) over a domain of size s can be approximated by: 


F(7 


a) 


/A (Apiasma'S/2) 

^ (Apiasn.aS/2)^ 

Tj I , sin (Apjasma'5/2) 

7o( 7 o)- - -^ 

(AplasmaS/2)' 


Notice then that the oscillation wavenumber Aqsc is given by 


A 


OSC 




A 


plasma? 


( 21 ) 


( 22 ) 


SO that the oscillation length is smaller than previously and is of the same order as 
the size of the domain of coherence considered. 
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Figure 1: Ratio between the probability of oscillation of a photon into a pseudoscalar 
including the effect of the intergalactic plasma and the probability of oscillation when 
this effect is not considered (as done in [^). The curves are drawn as a function 
of the photon energy, for Ma ~ 4 x lO^^GeV, nia ~ 10“^®eV and Bq ~ 10“®G, and 
for a distance of flight 0.5 Mpc (dashed line), 1 Mpc (solid line) and 2 Mpc (dotted 
line). 


It follows that the probability of oscillation is lower than in the previous case 
(with no plasma effects taken into account) and that it is no longer achromatic (see 
Fig. §. Supernovae observations not only argue for a dimming of distant supernovae 
but also argue for an effect that is achromatic, we will discuss this in more details 
in the next section. 

It is important to realize that when plasma effects are considered there are two 
sources for the loss of coherence, spatial fluctuations in the magnetic held and varia¬ 
tions in the number density of electrons (i.e. changes in the plasma frequency). For 
example, in the case of photon-pseudoscalar conversion in the interstellar medium 
of our galaxy |T^, the coherence length is most likely set by the huctuations in the 
plasma frequency. When this is the source of coherence loss one expects generically 
that s in equation (pT]) also depends on frequency. 
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4 Chromaticity constraints and coherence length 


As we mentioned in the previous section, once plasma effects are considered the 
conversion probability can depend on photon frequency. In this section we consider 
such constraints. 

Supernovae observations put a constraint on what is called the color excess be¬ 
tween the B and V wavelength bands (-E[B — V]). The color excess is defined 
as 


E[B - 1/] = -2.5 


logio 


F°{B) F%Vy 
F^B) F°{V)\ ’ 


(23) 


where F° (resp. F®) is the observed (resp. emitted) flux and the B (resp. V) band 
corresponds to 0.44pm (resp. 0.55pm). Observations constrain F[B — V] to be lower 
than 0.03 |^. This can be translated to 


F (7 ^ a)v 


F(7 ^ a)B 
F (7 ^ a)v 


< 0.03, 


(24) 


or equivalently to the statement that the conversion probability has to scale with 
photon wavelength weaker than near the visual band. 

We now consider different limits of equation (|^) to investigate the chromatic 
behavior. It is useful to introduce the plasma length scale £p = We can 

rewrite (|^) as 


F (7 ^ a) ~ {2AMipysm^{s/2ip). (25) 

In any astrophysically realistic situation, the effective coherence length is either 
set by the magnetic held or by the plasma frequency, so that it depends on proper¬ 
ties of the IGM which are expected to have a signihcant dispersion^, and will then 
similarly exhibit such a dispersion. 


We start by considering the case in which the coherence length is set by the 
magnetic held domains. For simplicity we will take s to be distributed as a Gaussian 
variable with mean s* and dispersion j3s*. The average of equation (|^) over s yields. 


(F(7 ^ a)) ~ 2 (Am4)' 


cos s 


74) e 


-(/3s*/G)V2' 


(26) 


In the limit s* S> ip, it reduces to 

(F(7 ^ a)) ~ 2 (Am4)'- (27) 

^For example, in the model of Ref. the magnetic field in the intergalactic medium was 
generated in quasars and expelled in their outflows. The distribution of bubble sizes at redshifts 
around zero is predicted to be extremely broad with two main components, sizes ranging from 
0.5 Mpc to about 5 Mpc, and a mean size of order 1 Mpc. The biggest bubbles have the largest 
magnetic field. 















The plasma length ip scales with the photon frequency so that the conversion prob¬ 
ability is proportional to which is ruled out by the constraint in equation (^). 
Notice also that a conversion probability that grows with frequency would produce a 
very significant reddenning of sources located at distances such that the probability 
reaches its saturation value for the shorter wavelengths. The absence of such red¬ 
denning may provide even more stringent constraints on the required achromaticity 
than that of Eq. (24). 

In the opposite limit, s* -C ip, we get 

(F(7 ^ a)) ~ {Ams*Y, (28) 


which is achromatic if s* does not depend on frequency. The total probability of 
conversion will still be given by Eq. with s* replacing s, however since one 
must have s* -C ip for (^8]) to hold, and since ip ~ 0.1 Mpc, this means that the 
total probability of conversion will be lower than the one computed in [|^ by at 
least a factor 10 to 100. One can try to overcome this by changing the coupling 
of the pseudoscalar, it is however difficult because the coupling is already near the 
astrophysical bound. Alternatively, the conversion probability may remain achro¬ 
matic as well as sufficiently large if the intergalactic magnetic field were stronger in 
domains of size s* <C ip. Faraday rotation measurements of distant quasars impose 
a conservative bound Bq < 10“® G ^Mpc on the strength of an intergalactic 
magnetic field coherent over scales Sc (it may be somewhat stronger depending 
on its spatial structure ||^). The conversion probability could be achromatic in the 


visible band and become of order unity at cosmological distances if the intergalactic 
magnetic field had some very definite spectral properties, for instance if Bq were of 
order 10“®G over domains of average size of the order of 10 kpc, and sufficiently 
weaker on larger domains. Additional constraints on the spatial distribution of the 
intergalactic magnetic field compatible with the proposed mechanism arise from pre¬ 
venting excessive dispersion in the observed peak luminosity of distant supernovae. 
Glusters, for instance, may have magnetic fields significantly strong and extended 
to make the conversion probability of order unity as photons get across them. 

Another possibility is to consider the case s* ~ ip and to require that the average 
in equation ( [2B|) be achromatic enough. This usually does not happen, although it 
can be accommodated by correctly choosing (3 and s*, or by a tight correlation be¬ 
tween the strength of the magnetic field and the sizes of the domains, and requires 
very special, and likely unrealistic, statistical properties of the IGM 


Let us now turn to the case where the coherence length is determined by the 
spatial variations of the electron density. In this case the required achromaticity 
would even be more of a fine tuning as the coherence length s will also depend on 
frequency. Just as in the case of the interstellar medium studied in |^, the frequency 
dependence will be set by the clustering properties of the electron density. 

In the case at hand we expect the intergalactic medium to be very clumpy and 
complicated. Gurrent numerical simulations indicate that the baryons today can be 
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found in several phases. About 30 % by mass is in a warm phase (T ~ SOOOA') 
that hlls most of the volume of the universe, about another 30% is in a warm-hot 
phase that resides in non-virialized objects such a hlaments and the rest resides in 
virialized objects such as clusters and in condensed forms such as stars and cool 
galactic gas (see e.g. [Q and references therein). 

The most relevant phase for our study is the warm one because it hlls most of the 
volume. It could be clumpy on scales smaller than ip, so the dumpiness of the IGM 
may be the most likely source of coherence loss. However a detailed study of the 
loss of coherence should probably involve studying lines of sight across these type 
of simulations. It is important to realize that even in the limit where the clumping 
scale, Lc, of the warm phase of the IGM is smaller than ip we still expect that s 
will depend on frequency. For coherence to be lost, the random component of the 
accumulated phase of the oscillation, 0 , has to be of order one. The phase on each 
segment of length Lc is L^/ip and accumulates over different segments as a random 
walk, 0 ~ 1. In this limit, we estimate s ~ i^pjLc which will again 

induce a chromaticity that is ruled out by observations. 


It seems that the only natural ways to avoid the chromaticity constraint is (i) 
to assume that the magnetic held is responsible for setting the coherence length 
of the oscillation and assume that s* <C ip, in which case either the coupling of 
the pseudoscalar needed to accommodate the dimming of the supernovae becomes 
uncomfortably large or the magnetic held must have very dehnite strength and 
spectral features, or (ii) to have very constrained properties of the IGM. 

Finally, we note that the situation is not improved by giving a higher mass to the 
pseudoscalar in order to have > Apiasma; this choice leads as well to chromaticity 
(since and Apiasma have the same spectral dependence) and lowers the probability 
of oscillation (and the mixing angle). A last logical possibility is that the mass of the 
pseudoscalar is such that A^ and Apiasma are of the same order. In this case, a strong 
mixing regime is possible whenever the electron density is such that A^ and Apiasma 
coincide with each other with accuracy Am- This can happen in the IGM from 
the statistical fluctuations in the electronic density. When this is the case, a strong 
photon-pseudoscalar conversion can in principle take place (this is very analogous 
to the resonant MSW effect of neutrino physics). However, for the transition to be 
signihcant one has to maintain the resonant condition (A^ ~ Apiasma with accuracy 
Am) over a distance of the order of the oscillation length, tt/Am, which given the 
numerical values ([T^) is also very unlikely (not to mention the fact that, a dimming 
of SNIa induced by such a resonant conversion would lead to a large dispersion in 
the observed SNIa magnitude). 
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5 Conclusions 


We have shown that one can not ignore the effect of the intergalactic plasma to 
derive how the Inminosity of distant sonrces, snch as snpernovae, is affected by a 
mixing with a hypothetical psendoscalar particle. In most of the parameter space, 
this effect either renders the oscillation frequency-dependent or lowers too much the 
oscillation probability. There is a slight hope to accommodate the mechanism of 
T5| if the IGM has very specihc statistical properties. 
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